Exact solution of mean geodesic distance for Vicsek fractals 
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The Vicsek fractals are one of the most interesting classes of fractals and the study of their 
structural properties is important. In this paper, the exact formula for the mean geodesic distance of 
Vicsek fractals is found. The quantity is computed precisely through the recurrence relations derived 
from the self-similar structure of the fractals considered. The obtained exact solution exhibits that 
the mean geodesic distance approximately increases as an exponential function of the number of 
nodes, with the exponent equal to the reciprocal of the fractal dimension. The closed-form solution 
is confirmed by extensive numerical calculations. 
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The concept of fractals plays an important role in 
characterizing the features of complex systems in nature, 
since many objects in the real world can be modeled by 
fractals In the last two decades, a great deal of activ- 
ity has been concentrated on the studies of fractals @, Q . 
It has been shown @, i, @, 0, IE H, Efll, E3 that regular 
fractals capture important aspects of critical percolation 
clusters, aerogels, amorphous solids, and unusual phase 
transition in the Ising model. Among various regular 
fractals, the Vicsek fractals [l^] are a class of typical can- 
didates for exact mathematical ones and have received 
much attention. A variety of structural and dynamical 
properties of Vicsek fractals have been inve stig ated in 
much detail, including eigenvalue spectrum [13j |. eigen- 
states 11411. L aplacian spectrum [15| . random walks [16j . 
diffusion jl7j ], and so on. The results of these investi- 
gations uncovered many unusual and exotic features of 
Vicsek fractals. 

A central issue in the study of complex systems is 
to understand how their dynamical behaviors are influ- 
enced by underlying geometrical and topological prop- 
erties [l8l Il9j . Among many fundamental structural 
characteristics {20|, mean geodesic distance is an impor- 
tant topological feature of complex systems that are of- 
ten described by graphs (or networks) where nodes (ver- 
tices) represent the component units of systems and links 
(edges) stand for the interactions between them [2l|, • 
Mean geodesic distance is defined as the mean length of 
the shortest paths between all pairs of nodes. It has 
been well established that mean geodesic distance di- 
rectly relates to many aspects of real systems, such as 
signal integrity in communication networks, the prop- 
agation of beliefs in social networks or of technology 
in industrial networks. Recent studies indicated that a 
number of other dynamical processes are also relevant to 



mean geodesic distance, including disease spreading [23| , 
random walks [24j |. navigation [25| . to name but a few. 
Thus far great efforts have been made to valuate and 
understand the mea n g eodesic distance of different sys- 

tems m, m m m m mi . 

Despite the importance of this structural property, to 
the best of our knowledge, the rigorous computation for 
the mean geodesic distance of Vicsek fractals has not 
been addressed. To fill this gap, in this present paper 
we investigate this interesting quantity analytically. We 
derive an exact formula for the mean geodesic distance 
characterizing the Vicsek fractals. The analytic method 
is on the basis of an algebraic iterative procedure ob- 
tained from the self-similar structure of Vicsek fractals. 
The obtained precise result shows that the mean geodesic 
distance exponentially with the number of nodes. Our 
research opens the way to theoretically study the mean 
geodesic distance of regular fractals and deterministic 
networks [ID, )H, [Hj] . In particularly, our exact solution 
gives insight different from that afforded by the approx- 
imate solution of stochastic fractals. 

The classical Vicsek fractals are constructed itera- 
tively [II El. We denote by V flt (t > 0, / > 2) the Vic- 
sek fractals after t generations. The construction starts 
from (t = 0) a star-like cluster consist of / + 1 nodes ar- 
ranged in a cross- wise pattern, where / peripheral nodes 
are connected to a central node. This corresponds to Vf t o- 
For t > 1, Vf,t is obtained from Vf t t—\- To obtain Vf t i, 
we generate / replicas of V/,o and arrange them around 
the periphery of the original V/.o, then we connect the 
central structure by / additional links to the corner copy 
structures. These replication and connection steps are re- 
peated t times, with the needed Vicsek fractals obtained 
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whose fractal dimension is 
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Fig. [TJ we show schematically the structure of V^- Ac- 
cording to the construction algorithm, at each time step 
the number of nodes in the systems increase by a factor 
of / + 1, thus, we can easily know that the total number 
of nodes (network order) of V/,t is N t = (/ + l) t+1 . 
After introducing the Vicsek fractals, we now investi- 



FIG. 1: Illustration of a particular Vicsek fractal V-3,2- 



gate analytically the mean geodesic distance between all 
the node pairs in the fractals. We represent all the short- 
est path lengths of Vft as a matrix in which the entry 
dij is the geodesic distance from node i to node j, where 
geodesic distance is the path connecting two nodes with 
minimum length. The maximum value D t of dy is called 
the diameter of V/^. A measure of the typical separation 
between two nodes in Vf t t is given by the mean geodesic 
distance L t defined as the mean of geodesic lengths over 
all couples of nodes: 

t = St m 

1 N t (N t - l)/2 ' 1 ' 

where 

St= £ dy (2) 

denotes the sum of the geodesic distances between two 
nodes over all pairs. 

We continue by exhibiting the procedure of the deter- 
mination of the total distance and present the recurrence 
formula, which allows us to obtain St+i of the t + 1 gen- 
eration from St of the t generation. By construction, the 
fractal Vf^+i is obtained by the juxtaposition of / + 1 

copies of Vf.t that are consecutively labeled as vj 1 ^, Vj I 2 ) , 

• • • , V^ t +1 \ see Fig. [2J This obvious self-similar struc- 
ture allows us to calculate St analytically. It is easy to 
see that the total distance St+i satisfies the recursion 
relation 

S t+1 = (f + l)S t + Q t , (3) 

where 0t is the sum over all shortest path length whose 
endpoints are not in the same V/j branch. The solution 
of Eq. © is 

t-i 

S t = (f+ I f So+Y, [(/ + lY^-'Qrn] ■ (4) 

m=0 



FIG. 2: A schematic illustration of the iterative construction 
for Vf,t+i, which is obtained by joining / + 1 copies of Vf t t 
denoted as Vj 1 ^ , V?^ , ■ ■ ■ , Vj^, and Vj t , respectively. 



Thus, all that is left to obtain 5* is to compute O m . 

The paths that contribute to Ot must all go through 
at least one of the 2/ edge nodes (such as G, X, Y, 
and Z in Fig. \2§ at which the different Vf t t branches are 
connected. The analytical expression for t , named the 
crossing path length, can be derived as below. 

Denote O"'' 3 as the sum of all shortest paths with end- 
points in Vj°^ and vj^. For convenience, we denote by 

Vf ^ the central branch of Vf t t+i. According to whether 
or not the two branches are adjacent, we sort the cross- 
ing path length Qf'@ into two classes: Q 1 ^ (0 > 1), Qf' 6 
(cp > 1, 6 > 1, and ip ^ 9). For any two crossing paths 
in the same class, they have identical length. Therefore, 
in the following computation of @t, we will only consider 
O]' 2 and Oj' 3 . The total sum O f is then given by 

e t = /xe t 1 - 2 +Q xef, (5) 

To calculate the crossing path length O]' 2 and Q 2,3 , 
we give the following definition and notations. We define 
external nodes of V/,< as the nodes that will be linked to 
one of its copes at step t+1 to form Vf+i t t- Let dt denote 
the sum of length of the path from an external node of 
Vf t to all nodes in Vfj including the external node itself. 

We assume that the two branches vj^ and V* 2 ) are con- 
nected at two nodes X and G, which separately belong 
to V$ and V$, and that V$ and V$ arc linked to 

each other at two nodes Y and Z that are in vj^ and 

Vjj , respectively. 

In order to determine dt, we should compute the di- 
ameter D t of Vf t t first. By construction, one can see that 
the diameter D t equals the path length between arbi- 
trary pair of external nodes of Vfj- Thus, we have the 
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following recursive relation: 

A+i = 3A+2. 



(6) 



Considering the initial condition Do = 2, Eq. |6j) is solved 
inductively to obtain 



A = 3 



t+i 



1. 



(7) 



which is independent of /. 

We now calculate the quantity d t +\. Let K denote the 

(2) 

external node of Vf t t+i, which is in the branch Vi t . By 
definition, dt+i can be given by the sum 

dt+i = Y^ dxj 

iev>, t+ i 

= Y dKu+ E djfu + (/-!) Y d 



Kw 



uev 



w€V t 



(3) 



Y d Ku + (f-i) y d 



/%»■■ (8) 



toev. 



We denote the second and third terms in Eq. {SJ by <?t 
and q t , respectively. Thus, d t+ i = d t + g t + q t . The 
quantity g t is evaluated as follows: 



(i) 



- d t + iV* x (A + l), (9) 
where dux — A and dox = 1 were used. Analogously, 
qt = (/-I) ( d ^G + ^gx + dxy + rfyz + rfz lu ) 

= (/-l)[dt+JV t x2(A + l)]- (10) 
With Eqs. © and HQ]), Eq. © becomes 

dtf l = (/ + l)dt + (2/ - 1) x N t x (A + 1). (11) 

Using JV t = (/ + 1)*+!, A = 3 t+1 - 1 and do = 2/ - 1, 
Eq. (|11|) is resolved by induction 



rf, = -(2/-l)(3 t+1 -l) (1 + /)*. 



(12) 



With above obtained results, we can determine the 
length of crossing paths ®]' 2 and 6 2 ' 3 , which can be ex- 
pressed in terms of the previously explicitly determined 
quantities. By definition, O]' 2 is given by the sum 

Y di 3 



(ckx + dxG + d G j) 
N t Y d iX + (N t ) 2 + N t Y d Gj 
2^V t 5] doc + (JVt) a , (13) 



where we have used the equivalence relation 
Proceeding similarly, 



E 



^ (d iG + d G x + 4r + rfyz + dzj) 
27V t J] d lG + 2(7V t ) 2 



(iV t ) 2 dxy 



i€V, 



(2) 



27V t ^ d lG + (7V t ) 2 (A + 2). 



(14) 



iev 



Inserting Eqs. (fT3]) and (fl4|) into Eq. (El, we have 



e t = (f+f)N t dt+f (N t y + 



2 , /(/"I) 



(-/V t ) 2 A- (15) 



Substituting Eq. (fl~5|) into Eq. ([4]) and using the initial 
value So = f 2 , we can obtain the exact expression for the 
total distance 



St 



1 



(/+ir / 2 (/+i)*(3 t+i +i) 



6/ + 4' 

3(3 t+1 -l)/ 8 (/ + l) t +4((/ + l)*-l) 
2/ (3 1+i (/ + l)*- 4(1 + /)* + !) 



(16) 



Then the analytic expression for mean geodesic distance 
can be obtained as 



Li 



1 



(3/a + 5/ + 2)[(/ + l)*-l] 
+ 3 (3 t+1 - 1) f(f + if + 4 ((/ + 1)' - 1) 

- 2/ (3 1+t (/ + l) f - 4(1 + /)* + !) 



/ 2 (/ + l)'(3 w + l) 
(17) 



In the infinite system size, i.e., t — > oo 
L t ~3 t+1 = (A^WTiy, 



(18) 



where the exponent fcnqzxl is equal to the reciprocal of 
the fractal dimension. Thus, the mean geodesic distance 
grows exponentially with increasing size of the system. In 
contrast to many recently studied network models mim- 
icking reahlife systems in nature and society [U, [12] , the 
Vicsek fractals are not small worlds despite of the fact 
that these fractals show similarity (fractality) observed 
in many real-world systems. 

We have checked our analytic result against numerical 
calculations for different / and various t. In all the cases 
we obtain a complete agreement between our theoretical 
formula and the results of numerical investigation, see 
Fig. El 

To sum up, in complex systems the mean geodesic dis- 
tance plays an important role. It has a profound impact 
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FIG. 3: Mean geodesic distance Lt versus network order Nt 
on a log- log scale. The solid lines are guides to the eyes. 

on a variety of crucial fields, such as information pro- 
cessing, disease or rumor transmission, network design- 
ing and optimization. In this paper, we have derived an- 



alytically the solution for the mean geodesic distance of 
Vicsek fractals which have been attracting much research 
interest. We found that in the infinite network size limit 
the mean geodesic distance scales exponentially with the 
number of nodes. Our analytical technique could guide 
and shed light on related studies for deterministic frac- 
tals and network models by providing a paradigm for 
calculating the mean geodesic distance. Moreover, as a 
guide to and a test of approximate methods, we believe 
our vigorous solution can prompt the studies on random 
fractals. 
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